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Abstract 

The QED processes e+e~ —>■ e'^e'fi^fi" and e+e^ e+e^r+r^ are studied with 
the L3 detector at LEP using an untagged data sample collected at centre-of-mass 
energies 161 GeV < ^/s < 209 GeV. The r-pairs are observed through the 
associated decay of one r into euu and the other into ttttz/. The cross sections are 
measured as a function of y/s. For muon pairs, the cross section of the 77 — > ^~ 
process is also measured as a function of the two-photon centre-of-mass energy for 
3 GeV < PV^^ < 40 GeV. Good agreement is found between these measurements 
and the 0{a'^) QED expectations. In addition, limits on the anomalous magnetic 
and electric dipole moments of the r lepton are extracted. 



Submitted to Phys. Lett. B 



1 Introduction 



The pair production of charged leptons in two-photon colhsions offers a unique opportunity 
to test QED to C(a^) over a wide kinematical range. The e+e~ — > e~ ji^ ji" and e+e~ — > 
e+e~r+T~ reactions are studied with the L3 detector [1] for untagged events, in which the e+ 
and e~, scattered at small angles, are not required to be observed. 

Figure 1 shows the lowest order processes which contribute to this final state: multipe- 
ripheral, bremsstrahlung, annihilation and conversion, for a total of 12 possible diagrams. For 
untagged events the multipcriphcral process dominates the cross section. The contribution of 
other processes is below 1%. 

The e"'"e~ — > e^e~t^l~ reactions, where £ — e, /j, or t, were previously studied for untagged 
two-photon events at e+e~ centre-of-mass energy, ^/s, close to the Z mass [2] . Good agreement 
was found between the measurements and the QED expectations. In this Letter, the production 
of yU-pairs is studied in the range 161 GeV < < 209 GeV and the production of r-pairs in 
the range 189 GeV < ^/s < 209 GeV. The individual energies and luminosities are listed in 
Table 1. As the e''"e~ — > q" jj/^ jjT channel benefits from high statistics, the cross section 
of the process 77 — ^ fi~^fi~ is also measured as a function of the two-photon centre-of-mass 
energy, VKy^. The e+e~ — > e+e~r+r~ analysis is restricted to the exclusive final state where 
T" — > e'ly^-Ug and t+ — > tt+tt^i/^-'^^ which arises from 9.07 ± 0.01% of all r-pair decays [3]. 

The process e+e~ — > e+e~T+T~ is also used to constrain the anomalous magnetic and electric 
dipole moments of the r lepton, as proposed in Reference 4. 

2 Data and Monte Carlo Samples 

The events are mainly accepted by the charged-particle [6] and the inner track triggers [5]. 
The former requires at least two charged particles with a transverse momentum pt > 150 MeV, 
back-to-back within an acoplanarity angle of ±41°. The latter is based on a neural network, 
has no requirement on the acoplanarity angle of the tracks and extends the acceptance from 
the polar region 30° < 9 < 150° to 15° < 6 < 165°. A fraction of the e^e^ — > c+c^/i+^~ events 
is also accepted by the muon trigger and a fraction of the e"'"e~ e'^e~T~^T~ events by the 
calorimetric energy trigger [7]. 

The DIAG36 [8] generator is used to calculate at 0{a'^) the full set of diagrams shown in 
Figure 1. To obtain the efficiencies of the e+e" — > e^e'jJL^jjr channel, high statistics samples are 
generated in the range 3 GeV < W^^ < 40 GeV, for each value of ^/s. The e+e~ e"'"e~r+r~ 
events are generated in the full phase space with the Vcrmaseren Monte Carlo [9], which takes 
into account only the dominating multiperipheral diagrams, shown in Figure la. 

For background studies, the following event generators are used: KORALZ [10] for the 
e+e- T+T- and e+e" fx+fj,' processes and LEPWW [11] and PYTHIA [12] for W and Z 
boson pair-production and decays into leptons, respectively. In the tau-pair analysis, exclusive 
hadronic two-photon processes are generated with EGPC [13] and inclusive hadron production 
with PHOJET [14]. 

All generated events are processed through the full L3 detector simulation based on the 
GEANT [15] and GHEISHA [16] programs and are reconstructed following the same procedure 
as for the data. Time dependent detector inefficiencies, as monitored during the data taking, 
are also included. 

Charge conjugate processes are included throughout this Letter. 
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3 Event Selection 



3.1 e"'"e e+e fi'^fi 

The muon pairs are selected using information from the central tracking chamber (TEC) and 
the muon spectrometer. The selection requires: 

• exactly two tracks with at least 12 hits each and opposite charges, having a distance of 
closest approach to the nominal interaction vertex in the plane transverse to the beam 
direction smaller than 5 mm; 

• two well reconstructed muons in the muon chambers corresponding to the charged tracks; 

• a fiducial volume | cos 6*^1 < 0.8, where 9^ is the angle between the muon and the beam 
axis; 

• the momentum of the muons between 2.5 GeVand 40 GeV; 

• muon tracks pointing to the primary vertex, with time-of-flight consistent with the beam 
crossing, in order to suppress background from cosmic rays, hadrons decaying in flight 
and punch-through hadrons; 

• a di-muon effective mass, M^^, which measures W-y-y, between 3 GeV and 40 GeV. 

The numbers of events selected at different ^/s are shown in Table 1 together with the 
selection and trigger efficiencies. The total background contribution, estimated by Monte Carlo, 
is below 1%, and consists mainly of events from the e+e~ — > e+e~T+T~, e''"e~ — > t^t~ and 
e'^e" — > processes and cosmic rays. The distributions of the di-muon effective mass and 

of the momentum of the higher energy muon are presented in Figure 2a and 2b together with 
the Monte Carlo predictions. The expected distributions agree well with the data. 

3.2 e+e- e+e-r+r" 

The selection of tau-pairs, through the associated decays r~ — > e~z/^z7p and r+ tt+tt^z/^, is 
based on information from the TEC and the electromagnetic calorimeter (ECAL). It requires: 

• a total energy in the calorimeters less than 40 GeV, to exclude e+e~ — > r+r" events; 

• exactly two charged tracks with at least 12 hits each and opposite charges, having a 
transverse momentum greater than 0.3 GeV, a distance of closest approach to the nominal 
interaction vertex in the plane transverse to the beam direction smaller than 10 mm and 
a corresponding ECAL signal; 

• two photons, defined as isolated showers in the ECAL with energy greater than 100 
MeV distributed over at least two crystals. There must be no track within 150 mrad 
around the shower direction and the ratio between the energies deposited in the hadronic 
and electromagnetic calorimeters must be less than 0.2. 
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The electron identification for the reaction t~ — > e~i/^z7g is based on an ECAL cluster, with 
a shower shape consistent with that of an electromagnetic particle, matching with a charged 
track within lOOmrad in the plane transverse to the beam direction. The momentum of the 
electron candidate must be greater than 600 MeV. To achieve high efficiency and high purity, 
the electron identification is based on a neural network [17] which combines ten variables: the 
energy in ECAL, the momentum, the ionisation energy loss in TEC, the ratio of the transverse 
energy in ECAL to the transverse momentum in TEC, the number of crystals in the shower, 
three inputs describing the shower shape in ECAL, the corresponding energy in the hadronic 
calorimeter and its fraction within a 7° cone. The electron identification with the neural network 
has an efficiency of 87.7 ± 0.2% with a purity of 94.7 ± 0.2%, as determined from Monte Carlo 
events. 

To identify r"*" tt+tt'^z/^ decays, we require the two photons to be compatible with a 7r°. 
The remaining charged particle is considered to be the 7r+ candidate. No additional selection 
cut is applied on the 7r+. The two-photon effective mass distribution in Figure 3a shows the 7r° 
peak. A gaussian fit to this peak gives a mass of 134.6 ± 0.6 MeV and a width of 6.8 ± 0.7 MeV, 
compatible with the expected detector resolution. We require the two-photon effective mass to 
be within the range from 115 MeV to 155 MeV. To reject exclusive final states, as for example 
e"'"e~ — > e"'"e~a2(1320) e+e~7r"'"7r~7r°, we require the total transverse momentum imbalance 
I to be greater than 0.2 GeV. Figure 3b compares the | distribution of data and 

Monte Carlo. The excess of data for | < 0.2 GeV is due to exclusive two-photon processes 

not included in the Monte Carlo. 

With these criteria, 266 events are selected. As expected for the r"*" — > n'^'K^jJ^. decay 
channel, the tt+tt*^ effective mass is consistent with the p meson mass, as shown in Figure 3c. 
The energy distribution of the electron candidate is shown in Figure 3d. All data distributions 
are in good agreement with Monte Carlo simulations. 

Table 1 shows the number of observed events together with selection and trigger efficiencies. 
The latter are evaluated directly from the data [18]. In the analysis, a two dimensional trigger 
efficiency correction, based on the highest momentum track and the azimuthal opening angle 
between the two tracks, is applied to each event. The main background in the sample is 26% 
and is due to tau-pairs decaying to other final states, where leptons or pions are misidentified, 
or additional pions are not detected. The background from the e'^e~ — > t+t~ process, from 
resonant final states and from hadron production in two-photon collisions is less than 4%. The 
background from beam-gas and beam- wall interactions is found to be negligible. 

4 Results 

4.1 e"'"e~ e^e~ ii~ 

The cross-section of the process e+e" e^e~ p^ pT for 3 GeV < W-^-f < 40 GeV is measured 
for I cos^^^l < 0.8 and extrapolated to the full angular range. The results are given in Table 2 
for different values of 

For lower luminosities the systematic uncertainties are dominated by the uncertainty on the 
trigger efficiency, around 3%. At higher luminosities the main uncertainty of about 1.5% arises 
from the limited Monte Carlo statistics. The uncertainty due to the event selection is estimated 
by varying the selection criteria for the data samples with high integrated luminosity and is 
less than 1%. 

The cross section for the full angular range, presented in Figure 4, shows the expected slow 
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rise as a function of ^/s and is in good agreement with the QED prediction, as calculated by 
DIAG36 Monte Carlo. The cross section of the process 77 yU+^u" is derived by measuring 
the cross section of the e+e" e+e~yU+/i^ process in nine W^^ bins and scaling it by the 
two-photon luminosity function [19]. The values obtained at different ^/s are consistent within 
a given W^^ bin, as shown in Table 3 and Figure 5a. Combined results for the full data sample 
are listed in Table 3 and shown in Figure 5b together with the QED predictions. A good 
agreement is observed. 



4.2 e+e — > e+e r+r 

The total r-pair production cross section is given in Table 2. The cross section is lower than 
the e"'"e~ — > e~^e~fi~^^~ cross section because of the r-pair mass threshold of 3.6 GcV. The 
main contributions to systematic uncertainties comes from the variation of the cuts on | Xl^l 
and the electron momentum, both between 4% and 5%. The total systematic uncertainty 
due to selection criteria is estimated to be between 7% and 9%. Other sources of systematic 
uncertainties are the determination of the trigger efficiency, the Monte Carlo statistics and 
the uncertainty on the background level; their combined contribution is below 3%. Figure 4 
compares the measured cross section and the 0{a'^) QED calculation. A good agreement is 
found. 



4.3 Anomalous Couplings of the Tau Lepton 

Photon couplings to the tau lepton are in general due to its electric charge, the magnetic dipole 
moment and the electric dipole moment. They can be described by a matrix element in which 
the usual 7'' term is replaced by [20] : 



2mr 2mr 

where the form factors Fi{q^), ^2(5^) and F^{q^), functions of the four-momentum squared, q^, 
of the photon, are related to the tau charge, magnetic and electric dipole moments as: 

respectively. In the Standard Model, at tree level, Fi{q^) — 1 and F2{q^) — F^{q'^) — 0. 
Limits on -^2(0'^) and F^[q^) were derived from the decay width r(Z r+r^), relating the Zrr 
coupling to the photon couplings via SU{2) x U{1) invariance [21]. Direct studies of the 7rr 
couplings were performed at the Z pole, by the L3 [22] and OPAL [23] collaborations through 
the e+e~ — > Z — > t^t~^ process, and at the T(4S) by the BELLE collaboration through the 
e+e" ^ 7* — > T+T~ process [24] . 

Tau-pair production in two-photon collisions is sensitive to possible anomalous couplings of 
the tau lepton. Values of F2{q^) and dj. different from zero would modify the cross section of the 
e+e~ — >• e^e~T^T~ process [4]. By comparing the measured cross section with predictions [4] 
as a function of ^2(5^) and dr we obtain: 

|F2(0)| < 0.107, \dr\ < 1.14 • 10"^^ ecm 

at 95% confidence level, where the limit on each coupling is derived fixing the other coupling 
to zero. These bounds, limited by the size of the data sample, are in agreement with the more 
stringent published ones [22-24] and are derived from a different process. 
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/ a\t [pi)] 


^/ m 




No 


Nb 




161 


10.2 


18.4 ± 0.5 


99.4 ± 0.6 


193 


4 




172 


9.7 


18.9 ± 0.5 


98.4 ± 0.8 


223 


7 




183 


54.2 


18.4 ± 0.3 


99.7 ± 0.2 


1188 


15 




189 


170.3 


20.1 ± 0.3 


99.6 ± 0.1 


4025 


33 




196 


154.0 


18.9 ± 0.3 


99.7 ± 0.1 


3491 


36 




206 


192.7 


19.1 ± 0.2 


99.7 ± 0.1 


4576 


45 




189 


172.1 


1.18 ± 0.04 


71.8 ± 1.3 


85 


25 




196 


220.9 


1.29 ± 0.05 


60.1 ± 1.6 


97 


31 




206 


215.1 


1.08 ± 0.04 


58.0 ± 0.9 


84 


29 



Table 1: Centre-of-mass energies and corresponding integrated luminosities. The selection 
efficiency, e^, and trigger efficiency, etrig, are also given together with the number of observed 
events, Nd, and the background contribution, Nb- 





i^s) [ GeV] 


O-DATA [Pb] 


O-QED [pb] 


O-DATA [pb] 


o-QED [pb] 






cos6'^ < 0.8 


cos6'^ < 0.8 








161 


101.4 ± 7.2 ± 2.6 


115.4 


587 ± 43 ± 22 


668.3 




172 


119.2 ± 7.6 ± 3.1 


116.6 


700 ± 46 ± 27 


684.9 




183 


117.7 ± 3.4 ± 1.9 


118.3 


697 ± 20 ± 9 


700.7 




189 


117.1 ± 1.8 ± 1.8 


118.9 


697 ± 11 ± 9 


708.6 




196 


118.9 ± 2.0 ± 2.2 


120.3 


713 ± 12 ± 12 


717.8 




206 


122.6 ± 1.8 ± 1.7 


121.3 


738 ± 11 ± 8 


730.0 




189 






459 ± 68 ± 33 


442.6 




196 






454 ± 67 ± 42 


452.3 




206 






459 ± 76 ± 35 


466.0 



Table 2: The cross sections of the processes e+e"" e^e~ jjf^ iJi~ and e+e" e^e~T^T~ with 
their statistical and systematic uncertainties at different ^/s values compared to QED [8] ex- 
pectations. The cross section for e+e~ — * ii'^ ii~ for 3 GeV < W^^ < 40 GeV is given for 
both I cos 6*^1 < 0.8 and for the full solid angle. 
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w 

\GeV] 


ct(77 ^ ) [nb] 


\' — loo vjLV 






^UU VjrLV 


loo — zuy vjLV 




3- 4 


24.3 ±9.9 


28.0 ±6.6 


25.2 ±6.4 


27.7 ±6.1 


25.9 ±4.2 


26.8 


4- 5 


21.5 ±3.7 


23.0 ±2.7 


25.2 ±2.8 


24.9 ±3.0 


22.6 ± 1.9 


21.5 


5- 6 


18.4 ± 1.9 


18.6 ± 1.5 


21.6 ± 1.6 


19.1 ± 1.6 


18.7 ± 1.1 


18.6 


6- 7 


14.5 ± 1.5 


16.8 ± 1.3 


18.8 ± 1.4 


16.1 ± 1.3 


15.9 ±0.9 


17.0 


7- 8 


12.3 ± 1.5 


15.3 ± 1.3 


14.9 ± 1.4 


18.5 ± 1.7 


14.9 ± 1.0 


15.2 


8-10 


11.5 ± 1.3 


12.9 ± 1.0 


12.4± 1.1 


12.9 ± 1.1 


12.4 ±0.7 


13.2 


10 - 15 


8.9 ±1.0 


9.3 ±0.8 


9.1 ±0.8 


8.3 ±0.7 


8.9 ±0.5 


9.6 


15 - 20 


6.0 ±1.0 


6.1 ±0.7 


6.2 ±0.8 


6.6 ±0.8 


6.2 ±0.5 


6.2 


20-40 


3.1 ±0.6 


3.3 ±0.4 


3.2 ±0.5 


3.6 ±0.5 


3.3 ±0.3 


3.2 



Table 3: The cross section of the process 77 — * with its combined statistical and system- 

atic uncertainties as a function of W^^ for four different y/s values and their average together 
with the QED [8] expectations. 
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Figure 2: Distributions for selected e+e~ — > e'^e~ ii'^/j," events of a) the di-muon effective mass, 
Mfj,^, and b) the momentum of the most energetic muon, p^. The data are compared to the 
sum of the DIAG36 Monte Carlo and of the expected background, normalized to the integrated 
luminosity. 
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Figure 3: Distributions for r — > TTTTi/ candidates of a) the effective mass of final state photons, 
rrij^, b) the sum of the transverse momenta of the charged particles, | X^Ptl' ^) effective 
mass of the two pions, TOttttO- d) Distributions of the energy of the electron for r ei/i/ 
candidates. The data are compared to the sum of the Vermaseren Monte Carlo e+e~ — >■ 
e^e^r+r^ and of the background, normalized to the integrated luminosity. Arrows in a) and 
b) indicate the position of the cuts on the plotted variable, when all other selection cuts are 
fulfiUed 
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Figure 4: The cross section of the c+c^ - 
total cross section of the e"^e~ — > e+e^r+r^ process for W^^ > 2m,-. The data are compared 
to the QED calculations of DIAG36. The inner parts of the error bar represent the statistical 
uncertainties, the outer parts the systematic uncertainties 
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Figure 5: The cross section of the process 77 — > as a function of the 77 centre-of-mass 

energy for a) different values of ^/s and b) their combination. The data are compared to the 
QED calculations of DIAG36. 
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